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ABSTRACT 



An algorithm that extracts aerosol optical depth and particle size index from satellite 
observations of upwelling radiance over clear ocean areas was modified and studied. In 
order to examine the algorithm's performance on a regional scale, a previously analyzed 
data set, retrieved by the N’OAA-9 AVIIRR sensor, was reprocessed. The area of the 
study was in the central Pacific Ocean during the RITS-8S cruise from 7 April to 5 May 
1988. The results were compared to those of the earlier study and used to investigate 
changes caused by modifications to the algorithm. One significant modification was to 
account for the absorption of column water vapor in the calculations. In addition, the 
correlation between aerosol optical depth and the two-term Henyev-Greenstein phase 
function was investigated by normalizing measured radiance values. 

The results compared well with the earlier study and supported the modifications to 
the algorithm. Based on this study, aerosol optical depth was shown to be correlated to 
the Henyey-Greenstcin phase function. This knowledge helped in formulating improve- 
ments to the phase function. 

It was concluded that satellite imaging and processing on a regional scale is a useful 
way to study marine aerosols, that results can be improved by considering the effects of 
atmospheric water vapor absorption and other modifications to the calculations, and 
that there is room for improvement to the two-term Henyev-Greenstein phase function. 
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I. INTRODUCTION 



A. THE NATURE OF THE PROBLEM AND BACKGROUND 

As the world's population and knowledge have grown, so too has an interest in its 
climate. The ability to study our environment and anticipate possible climate changes, 
whether caused by humans or not, has far-reaching implications for future generations. 
Their very survival may depend not only on our ability to understand any changes, but 
also on what we do about them, whether it be to adapt to the changes or attempt to 
influence our climate. How well we might adapt to changes or whether we can influence 
our climate is still not known but we have begun to understand it by studying the overall 
radiation budget of the earth. 

Prior to the advent of the satellite, estimates of radiation entering and leaving the 
earth-atmosphere system were made using radiative transfer equations (Rao et al. 1990). 
Many factors had to be considered and incorporated into the model. Radiation in the 
atmosphere is attenuated by many processes. It is absorbed and scattered by molecules 
of air, aerosols, and meteorological features such as clouds, precipitation and fog. Aer- 
osols are solid and liquid particles, approximately one micron or less in diameter, sus- 
pended in air. This includes smoke, dust, haze, and some clouds. While the results of 
the radiation calculations were helpful, there are still many uncertainties. For example, 
in recent years there has been increasing interest in the warming of the earth's climate 
caused by increased concentrations of the so-called greenhouse gases, notably C0 2 , 
which trap outgoing long wave radiation. 

It was and still is necessary to parameterize the less well understood contributions 
of such features as aerosols and clouds. These features are inherently more difficult to 
measure than C0 2 because they have shorter time scales and they arise on a local or re- 
gional scale leading to higher variability. Sources of aerosols include natural, like 
volcanic aerosols and desert dust, anthropogenic, or human made, and biological, where 
gaseous metabolic byproducts are precursors to aerosols. McCormick and Ludwig 
(1967) proposed that a buildup of aerosols could have a cooling effect on the earth. 
Charlson and Pilat (1969) concluded that particle size distribution as well as scattering 
and absorption cross sections also played a role and could lead to warming in some 
cases, especially for anthropogenic aerosols. Coakley et al. (1983) studied naturally oc- 
curring tropospheric aerosols and found that they result in global surface cooling of 2-3 
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°C, an amount close to the increase in temperature due to such gases as CH 4 and X 2 0, 
or a doubling of C0 2 . In addition, the cooling effect of aerosols can be multiplied when 
they act as cloud condensation nuclei (CCN) resulting in the formation of cloud drop- 
lets. Clouds reflect incoming solar radiation leading to a cooling effect. Opposing this 
effect is a tendency for non-nucleated anthropogenic aerosols to absorb radiation 
thereby reducing cloud reflectance which leads to warming. Twomey (1977) and 
Twomey et al. (1984) concluded that for thin and medium optically thick clouds over 
continents, an increase in anthropogenic aerosols results in cooling. Thin and medium 
optically thick clouds are most common in the atmosphere. Thus, fluctuations in aer- 
osol concentration could play a role in the global climate with an increase in aerosol 
concentration contributing to a cooling effect. 

Furthermore, Charlson et al. (1987) has expanded on hypotheses of Shaw (1983) and 
later Lovelock (1986) that a biological feedback mechanism where marine 
phytoplankton, through production of marine aerosol precursors, play a role in climate 
regulation. 

Figure 1 is a schematic diagram of how the feedback mechanism works. In the di- 
agram the rectangles are measurable quantities and the ovals arc processes linking the 
rectangles. A positive sign in the oval indicates the effect a positive change in the pre- 
vious rectangle has on the following rectangle. A negative sign has a negative effect on 
the following rectangle for a positive change in the previous rectangle. 

Charlson et al. (1987) explain that the dimethylsulfide (DMS) produced by 
phytoplankton is oxidized in the atmosphere to form a non-sea-salt sulfate (XSS-SO^~). 
These XSS-SOj~'s act as CCX in the marine atmosphere. The argument is that DMS 
is produced in sufficient quantity to be the primary' source of CCX in the remote, un- 
polluted marine atmosphere. Since clouds form from CCX, and clouds act to reflect 
incoming solar radiation, phytoplankton might then play a role in regulating the earth's 
climate. The greatest uncertainty in the hypothesis is what effect any lowering of the 
earth's surface temperature and reduction of solar irradiance below the clouds would 
have on the phytoplankton and DMS production. 

Once satellites became available they were quickly tasked to help provide data on 
radiation or factors that influence radiation (Rao et al. 1990). For example, marine 
aerosols are atmospheric constituents that appear to play an important role in the 
overall radiation budget, as was shown. Hindman et al. (1984) describe experiments that 
establish a positive relationship between aerosol variations and satellite detected bright- 
ness variations in the marine boundary layer. In fact, direct measurements of marine 
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Figure 1. Biological feedback in the Earth's radiation budget (Charlson et al. 

1987): rectangles are measured quantities, ovals are linking processes. 

Plus or minus signs indicate the effect a positive change in a previous 
rectangle has on the following rectangle. 
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aerosol characteristics could be provided by ship or aircraft. While this information is 
valuable, it is expensive and difficult to obtain and only a small area of the world's 
oceans can be sampled and studied at a time. Satellites, on the other hand, can provide 
near continuous data over large areas. This capability is desirable for evaluating the 
effects of aerosols on a global scale. Meteorological satellites, however, have been de- 
signed to measure cloud albedo while sensitivities to radiances associated with aerosol 
optical depth variations are small (Durkee 1984). Thus, in order to determine aerosol 
characteristics, it is important that the satellite data be carefully processed. 

Aerosols can also play a role in military applications. Aerosols can scatter energy 
from a laser beam propagating through the atmosphere. Lasers can be used for satellite 
communications with submarines. The narrow beam of the laser is harder to detect than 
the normal radio communications now in use and is capable of higher data transmission 
rates to deeper depths. As submarines become more important as multimission plat- 
forms, it becomes more important to have secure communications with them (Painter 
19S9). In addition, line-of-sight communications operating in the near infrared can uti- 
lize aerosol scatter in the marine boundary layer to provide over-the-horizon ranges of 
30-300 miles (Mooradian 1981). 

Further investigations of aerosol characteristics arc clearly needed. The purpose of 
this thesis is to better utilize satellite data in order to more accurately measure aerosol 
characteristics so that their effect on the overall radiation budget of the earth and in 
other applications such as communications may be quantified. 

B. METHODOLOGY 

Durkee et al. (1991) describes an algorithm whereby multispectral radiance data 
obtained from the Advanced Very High Resolution Radiometer (AVHRR) onboard the 
U.S. National Oceanic and Atmospheric Administration (NOAA)-9 satellite is used to 
obtain estimates of aerosol optical depth and particle size index over clear ocean areas. 
By utilizing different satellite channel sensitivities, the particle size index can be meas- 
ured and the optical depth calculated. 

Modifications to the algorithm have been made. For example, Mahony (1991) 
studied and quantified the effects of column water vapor absorption on aerosol particle 
size index. This information can be used to improve the performance of the algorithm. 
In addition, a modification to account for variable earth-sun distance on solar irradiance 
was made. A change in the computer code to correct a mathematical error and an ad- 
justment to an over correction for ozone absorption were the final two modifications. 
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To study the modified algorithm on a regional scale, a data set used by Benedict 
(1989) was reprocessed. The results were compared to those of Benedict's (1989) earlier 
study. It was expected that all the modifications would improve the results with the 
water vapor correction being the most significant. An important extension of this pro- 
cedure was to study the correlation between optical depth and the two-term Hcnyey- 
Greenstein (HG) phase function used in the calculations. This was accomplished by 
studying a plot of normalized measured radiance versus scattering angle for various 
particle size indices. The purpose of doing this was to make optical depth independent 
of scattering angle; these variables might be related through the phase function. If there 
is a correlation, improvements to the phase function and therefore calculated optical 
depth are indicated. 

This methodology was chosen because the Benedict (1989) data set was available 
and provided a baseline from which to measure improvements. In addition, useful in- 
formation about the HG phase function had not yet been extracted from the Benedict 
(1989) data set. Once improvements to the algorithm are verified and permanently in- 
corporated into the computer code, new data sets can be studied with the goal of better 
understanding the role of DMS on the earth's climate. 

C. OBJECTIVES 

The specific objectives of this thesis are to compare the results to those of Benedict's 
(1989) study and note changes caused by modifications to the algorithm. In addition, 
improvements to the results derived from Mahony's (1991) correction for atmospheric 
water absorption will be studied. Finally, the form of the Henvey-Greenstein phase 
function obtained from this study is compared to that used by Frost (1988) in order to 
make recommendations for improvement. 
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II. THEORY 



A. RADIATIVE TRANSFER 

It will be shown in the following discussion that satellite measurements of upwelling 
radiance provide a means to determine aerosol optical depth and particle size informa- 
tion. These characteristics can help determine the amount and source of aerosols. The 
satellite wavelengths used in the approach are from the NOAA-9 AVHRR. Table 1 
gives channel bandwidths for reference. 



Table 1 . AVHRR CHANNEL BANDWIDTH: (in microns) from Kidwell (1986). 



Channel 


Radiance 


Bandwidth 


1 


Visible 


0.58 - 0.68 


2 


Red Near Infrared 


0.70 - 1.10 


3 


Thermal Infrared 


3.50 - 3.90 


4 


Thermal Infrared 


10.30 - 11.50 


5 


Thermal Infrared 


11.30 - 12.40 



In order to develop a practical form of the Radiative Transfer Equation (RTE) for 
use in marine aerosol calculations, it is necessary to review the theory and make some 
assumptions and modifications. From Liou (1980), the form of the RTE in a plane 
parallel atmosphere for scattering radiation is 

H = ^(<5. G) - ^ 'm G>(G, GW 

nF 0 p(Q-Q 0 )e~ 6l \ 



where 

L = Diffuse intensity of radiance. 

<5 = Optical depth. 

(i) 0 = Single scattering albedo. 

\x = cos 6 (observation zenith angle). 
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H„ = cos 6 0 (solar zenith angle). 

Q. = Solid angle ( 0 , (/>) (</> = azimuth angle). 

= Solid angle of direct radiation singly scattered into the path. 

Q' = Solid angle of multiply scattered radiation scattered into the path. 

p(£2, £2') = Scattering phase function from solid angle Q. to Q' (©). 

nF„ = Incoming radiative flux (solar flux density) adjusted for variable earth-sun dis- 
tance; different for channels one and two. 

The first term on the right represents the diffuse radiance attenuated from the beam 
by single scatter and absorption while the second and third terms represent contributions 
from multiple and single scatter into the beam, respectively. 

Optical depth is defined as the amount of absorbing material lying between two 
levels and can be represented as 



where o, xt is the extinction coefficient and H is the satellite height. The extinction coef- 
ficient, which is a sum of absorption and scattering contributions, is defined as 



where nr 2 is the particle cross sectional area, Q(m, )., r) is the extinction efficiency for 
cither scattering or absorption as a function of the complex index of refraction, the 
wavelength, and the particle radius, and dN{r)fdr is the slope of the aerosol size distrib- 
ution. 

For marine particles absorption is negligible when the wavelength is less than one 
micron. Then a ex „ which equals the sum of the extinction from scattering and absorp- 
tion, is nearly equal to o, C0 , (Shettle and Fenn 1979). The single scattering albedo is de- 
fined as o sca , divided by o tx , and is nearly equal to one. 

The scattering phase function, p(O), is the distribution of scattering by particles and 
is equal to one when summed over all solid angles. More will be said about the phase 
function below. 

If multiple scattering is neglected and only single scattering of direct solar radiation 
into the path is assumed, the RTE simplifies to (Liou 1980) 



*5 — f G ext dz 



ext Ud -’ 



( 2 . 1 ) 



o 




( 2 . 2 ) 
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L(<5; n, <f>) = L{5 X \ (i, <f>) exp[ - ((5i - S)jp] + 



:r ^ F oP (n, 0; - p 0 , 4 > 0 ) \\xp{ - [(<5' - d)ln + ^ I Pol) dS'lp, 

where 5, is the optical depth at the bottom of the atmosphere and S' is the increment 
of optical depth between the bottom of the atmosphere and the level of interest. 

Additionally, over ocean areas in the red and near infrared wavebands there is very 
little radiation emitted by the earth resulting in a further simplification of the RTE at the 
top of the atmosphere to (Liou 19S0) 

i(0; "• *) - TTTTV'’* 0 )! 1 - ex Pt - + IW1J. (2-3) 

For optically thin atmospheres, where optical depth is less than about 0.1, which is 
true over clear ocean areas, <5,(1/^ + 1 jp 0 ) 4 1 and equation 2.3 reduces to (Durkee 1984) 

Wp, 4>)* —^—p{Q)S v (2.4) 

Now the upwelling radiance can be broken down into its various contributions as 

L Total = L Molecules "b ^ Aerosols “b ^Sunglim ^ Ozone, 

where the subscripts stand for the contributions by molecules, or Rayleigh scatter, con- 
tributions from aerosols, contributions from sun glint, or specular reflection, and re- 
duction by ozone absorption, respectively. 

To get the aerosol radiance, the radiance measured by the satellite sensor was ad- 
justed by subtracting the Rayleigh scatter due to molecules following Turner (1973). 
Rayleigh scatter occurs when the size of the particle is much smaller than the wavelength 
of incident radiation and though strongly wavelength dependent, it does not vary spa- 
tially. 

An adjustment for sun glint is also necessary. Ramsey (1968) showed that sea sur- 
face reflectance is very small for red wavelengths (albedo = 0.5%) and zero for radiation 
with wavelengths greater than 0.7 microns, except where the sun-earth-satellite geometry 
results in sun glint. Therefore, geometries which could contribute to sun glint are 
avoided. Following Cox and Munk (1954), an angle of reflection is calculated from ge- 
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ometrv considerations. Then an elliptical area in the neighborhood of the reflection 
point where sun glint is possible is defined. This definition is based on the angular di- 
ameter of the sun. 

Since ozone absorbs weakly in channels one and two, a small ozone correction was 
added back. This correction was a modification to that of Benedict (1989). Now the 
aerosol radiance, L A , can be expressed as 

co F 

L a = 1(0; n, 4>) - L m - L sc + L 0 x -f^p(Q)S A , (2.5) 



where 5 A is aerosol optical depth. 

Griggs (1975) and Durkee et al. (1986) showed that there is indeed a linear re- 
lationship between aerosol radiance and optical depth. A method is now available to 
obtain aerosol optical depth from upwelling radiance by rearranging equation 2.5 to give 
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A ~ 



Mm 

MX©) ' 



(2.6) 



By taking the ratio of channel one to channel two aerosol radiance, given by 

l J Rpd (MXmM©)M) 

L\ir ( M ,\ir) 



a new parameter, the particle size index, is defined from which particle size information 
can be obtained: 



S 



12 “ 



Lr’ir 



1 $Red 
$ \IR 



(2.7) 



Since for smaller particles the aerosol dependent terms, a> 0i />(©), and <5 are larger at 
red wavelengths (Durkee 1984), smaller particles would cause S 12 to be larger. The in- 
crease in optical depth can be attributed to a steeper slope in the particle size distrib- 
ution as seen from equations 2.1 and 2.2. Thus the particle size index gives an indirect 
measure of the aerosol particle size distribution. For channels one and two the radiances 
are such that S 12 is normally greater than one. 

B. THE HENYEY-GREENSTEIN VARIABLE PHASE FUNCTION 

The scattering phase function describes the fraction of radiation scattered in a given 
direction. As discussed by Frost (1988), one single phase function is inadequate for 
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calculating optical depth when aerosols change over global scale areas. Therefore, a 
variable phase function, the two-term Henyey-Greenstein phase function (Lenoblc 
1985), was chosen and is given by 



p(Q) 



ftp -g?) 

(1 +g! 2 -2g, cos(0)) 



+ 



(1 - «)(1 -g 2 ) 

(1 + g\ + 2g 2 cos(6)) ’ 



( 2 . 8 ) 



where a is a weighting factor and g, and g 2 are asymmetry curve factors. All are empir- 
ical functions of S 12 (Frost 1988) and are given by 

a = 1.2- (0.25 xS 12 ), 



g i = 1.04 - (0.08333 xS 12 ), 



and 



g 2 = 1.2-(0.5S333 xS 12 ). 

Figure 2 shows a plot of the phase function versus scattering angle for various S 12 
values. Higher values of the particle size index (smaller particles) result in a curve that 
is closer to that of Rayleigh scattering. The MG phase function uses S 12 prior to cor- 
rection for Rayleigh scattering and ozone absorption. Since Rayleigh scattering, more 
in channel one, increases S 12 and ozone absorption, also more in channel one, decreases 
S, 2 . the difference is slight and the resulting particle size index and phase function do not 
vary greatly from corrected values (Durkee ct al. 1991). 

The HG phase function also varies with wavelength although it is only shown here 
for channel one. This introduces an error into the calculation of optical depth for 
channel two in equation 2.6. It will be shown in the analysis of results how this error 
may be corrected. 

C. WATER VAPOR CORRECTION TO THE PARTICLE SIZE INDEX 

The RTE in the above calculations assumed cloud free conditions. Another impor- 
tant consideration in atmospheric absorption is the moisture present in the atmosphere. 
Water vapor can produce as much as a 15% error in the phase function and aerosol 
particle optical depth (Mahony 1991). Dalu (1986) developed an algorithm to measure 
water vapor content using the channel four and five brightness temperature difference. 
The measured water vapor content of a column of air is used to determine the "dry" 
aerosol particle size index. 
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In Figure 3 the effect of water vapor can be seen in the transmittance (r) for 
AVHRR channels one and two. Note that for channel one t is greater than channel two 
indicating less absorption. Channel two radiance would tend to decrease more than 
channel one leading to overestimation of S I2 . This error could be misinterpreted as a 
change in aerosol properties or it could mask real changes. Thus, accounting for water 
vapor should improve the accuracy of the results. 

Mahony (1991) used the U.S. Air Force s Low Resolution Transmittance 
(LO\VTRAN)-7 atmospheric propagation model to study the effects of water vapor. 
Figure 4 shows the channel one and two relative response as well as the resulting radi- 
ances from Rayleigh scatter, Rayleigh and aerosol scatter, and Rayleigh plus aerosol 
scatter minus water vapor absorption. The radiances increase with smaller wavelength 
as the contribution from Rayleigh scatter increases. The radiance decreases with wave- 
length in both channels but more in channel two leading to overestimation of S 12 . 

An estimate of the water vapor content in an air column can be made using channels 
four and five of the AVHRR, at 11 and 12 microns respectively. Here absorption by 
water vapor is dominant while scattering is relatively insignificant. The benefits of using 
these channels are their collocation and temporal linkage with channels one and two, 
providing a water vapor snapshot along the same path as aerosol optical depth meas- 
urements. Although not the ideal sensor window for water vapor content, their use is 
justified within the accuracy of the approximations (Mahony 1991). 

The brightness temperature (T) is determined from the radiance received due to di- 
rect transfer from the ocean surface to the satellite sensor. Therefore the transmittance 
is directly proportional to the retrieved T. Only absorption reduces the radiance reach- 
ing the sensor, scattering is negligible. 

Dalu (1986) showed that by using the difference in brightness temperature between 
channels four and five the water vapor content could be obtained from 

w = A(T 4 — T 5 ) cos 0, 

where w is the water vapor content in grams per square meter and T 4 and T 3 are the 
brightness temperatures in degrees Kelvin for channels four and five, respectively, and 
6 is the scanning angle. The constant A is equal to 19600 grams per degree Kelvin per 
meter squared. 
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Figure 3. Water vapor effects on NOAA-7 channels one and two (Mahony 
1991): total transmittance is superimposed. 
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Figure 4. Water vapor influence on atmospheric radiation (Mahony 
1991): NOAA-7 response curves and LOWTRAN-7 used. 
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By comparing the water vapor content to the particle size index, Mahony (1991) 
obtained the following empirical relation for the particle size index corrected for water 
vapor: 



S, 2 (dry) 



£l_2 

(1 + . 0332 ^) ’ 



where w is the water vapor content in grains per square centimeter. This is the formula 
used in the calculations to obtain dry optical depth and particle size index. 



15 



III. PROCEDURE 



A. THE DATA SET AND THE SATELLITE 

The in-situ measurements were taken from the Research Vessel ( Rj V) 
Oceanographer, a NOAA funded research vessel, during the Radiatively Important Trace 
Species (RITS)-SS cruise conducted from 7 April to 5 May 19S8. The ship's track was 
in the central North and South Pacific Oceans from 50°N to 12°S along 170°W. Data 
measurements were obtained for every one degree of latitude and included DMS emis- 
sions in the sea and air, NSS-SOa - aerosol concentration, aerosol size distribution, and 
total particle concentration (Charlson and Bates 1989). Optical depth was also obtained 
using a hand held photometer. 

The satellite data were collected by the AVHRR sensor onboard the NOAA-9 sat- 
ellite. See Table 1 for AVHRR channel wavelength bands. The NOAA-9 is a sun- 
svnehronous, near polar-orbiting satellite with local equator crossing times of 0220, 
descending, and 1420, ascending. It has a nominal altitude of S33 km and resolution of 
1.1 x 1.1 km at sub- point. 

The satellite data utilized were for the same time frame and area as the ship's track 
from the digital archives of the National Environmental Satellite and Information Ser- 
vice (NESD1S). The 57 passes were stored as Global Area Coverage (GAC) data; GAC 
is an archiving technique that also reduces the resolution at sub-point to about 4.0 x 4.0 
km. 

B. PROCESSING FACILITIES 

The data were processed at the Interactive Digital Environmental (IDEA) Labora- 
tory at the Naval Postgraduate School in Monterey, California. The main computer 
was a Digital Equipment Corporation VAX 8250 with graphics display terminals. 
Processing included loading and reading the magnetic tapes to memory' and then proc- 
essing the data. The multi-channel analysis routine used to analyze the data was a 
FORTRAN program run on the VAX computer. Visual images were also produced in 
the IDEA lab and helped in the interpretation and analysis of the data. 

C. THE ALGORITHM 

The satellite data were analyzed using a scheme developed by Pfeil (1986) and 
modified by Frost (1988). In addition, corrections were made to solar irradiance to ac- 
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count for variations in earth-sun distance, for over correction of ozone absorption, and 
to correct the computer code for a mathematical error. Calculations were also made 
using the water vapor correction to S, 2 described by Mahony (1991). 

The pixels were subjected to a number of filters, described below, and composited 
in one degree by one degree boxes. To study the relationship between satellite and 
shipboard measurements, channel one aerosol optical depth and the particle size index 
were further filtered geographically to within five degrees longitude of 170° W and plus 
or minus three days of the ship's track. 

The first step in the algorithm is to test for valid data. After reading in the data and 
determining their geographical position, only data equatorward of 70°N or S is accepted. 
In addition, the data are required to be over ocean areas; here upwelling radiance from 
the surface is minimized. Next, sun angle is checked. Solar zenith angle must be less 
than 70° to avoid longer and more complicated optical paths. Data with specific ge- 
ometries in relation to the satellite are also eliminated to avoid sun glint (Cox and Munk 
1954). This eliminated a large portion of the western half of each pass. 

Before testing the pixels further, an adjustment to channels one and two were made 
to account for ozone absorption. This is of greatest concern for a wavelength of 0.69 
microns (Durkee et al. 1991). 

The calibrated pixels are then further processed to determine cloud contamination. 
If the brightness temperature of channel four is less than 273°K, the pixel is assumed to 
be contaminated by high clouds. Next, if the channel two albedo is greater than 40%, 
the pixel is assumed to be contaminated by low clouds. Then, if the channel one divided 
by channel two albedo ratio is greater than 1.5, the pixel is considered clear since cloud 
scattering is weakly dependent on wavelength. If the ratio is less than 1.5, the pixel 
undergoes further testing. The channel two albedo is compared to tnat of the pixels 
above, below, and to the left and right. If the difference between the maximum and 
minimum counts is less than five, the pixel is clear, otherwise it is partly cloudy. A count 
is approximately 0.1069% albedo (Kidwell 1986). 

Once a clear pixel has been determined, the Rayleigh radiance is subtracted from 

channels one and two (Turner 1973) and the optical depth and particle size index are 

calculated using equations 2.6 and 2.7. The derived parameters are: 

The Spatial Variables. Channel one and two radiances were used to obtain optical 
depth and particle size index. The "dry" values were also obtained using the algorithm 
of Mahony (1991). This information was composited temporally over the entire data 
set and was used to create images. 
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The Geometry Variables. These are variables such as normalized measured radiance 
as a function of scattering angle, latitude, and particle size index. They were com- 
posited geographically across all longitudes as well as temporally. One parameter, the 
normalized measured radiance as a function of scattering angle for different particle 
size indices, was used to study the two-term HG phase function. 

The Daily Variables. The particle size index and channel one optical depth both un- 
corrected and corrected for water vapor filtered to be within five degrees of longitude 
and plus or minus three days of the ship's track. This information was used for com- 
parisons with in-situ measurements. 

The calculated values were then stored in data files and used as needed to produce 
the necessary images and plots for comparison with those of Benedict (19S9) and the 
measurements of the R/V Oceanographer. 
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IV. RESULTS AND DISCUSSION 



The results are organized to show improvement to the analysis done by Benedict 
(1989), to show improvements to aerosol characteristics based on corrections for column 
water vapor as parameterized by Mahony (1991), and to investigate improvements to the 
HG variable phase function. 

Before looking at the results, it is useful to examine a summary of the synoptic sit- 
uation during the 7 April - 5 May 1988 time frame. The synoptic situation is important 
to the discussion since it helps identify sources of aerosol particles. Knowing the sources 
of the aerosols helps in interpreting the results. 

Representative Navy Operational Global Atmospheric Prediction System 
(NOGAPS) surface pressure and 500mb height patterns are shown in Figures 5 and 6. 
A generally westerly flow exists over the Asian continent in the 35-45°N latitude band 
indicating a source for natural and anthropogenic continental aerosols. Natural aerosols 
would tend to be larger in size than the marine aerosols being studied. Benedict (1989) 
discusses the occurrence of KOSA, or dust clouds from storms in the Gobi desert during 
April. Given typical translational rates (Shaw 1980), the dust from these storms could 
influence the aerosol concentration in the central Pacific. 

Another environmental factor discussed by Benedict (19S9) was the eruption of the 
Kilauea volcano in Hawaii (17°N, 157°W) producing a continual source of aerosols 
during this time frame. These aerosols were initially dominated by gaseous sulfur com- 
pounds which later changed to the solid phase in the atmosphere and grew as they mi- 
grated downstream. The debris from the eruption was carried across the path of the R/V 
Oceanographer's track by the Northeast Tradewinds and influenced the aerosol meas- 
urements. 

The composite optical depths for the entire experimental period are shown in Fig- 
ures 7 and 8. The optical depth is an indication of both the number and size of aerosols. 
The uncorrected optical depth is presented in Figure 7 which shows an area of relatively 
low optical depth in its lower right corner. A higher optical depth was expected due to 
the volcanic aerosols. The reason for the low optical depth is because most of the 
volcanic emissions are gaseous and very small and therefore do not greatly effect optical 
depth. It is also possible that the parameterization of the phase function is not correct 
leading to errors in optical depth. Figure 7 shows an area of high optical depth at the 
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Figure 5. NOGAPS 1200 UTC 10 April 1988 Surface Analysis (Benedict 1989) 
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Figure 6. NOGAPS 1200 UTC 10 April 1988 500n»b Analysis (Benedict 1989) 
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western edge of the area which extends eastward across the upper third of the figure. 
This indicates a potentially large number of particulates which could be either 
anthropogenic or natural in origin. When the particle size index is examined below, 
further information about the source of the aerosols may be inferred. From Figure 2 it 
is seen that for a given radiance, when the particle size index decreases the phase func- 
tion decreases. This leads to an increase in optical depth for most scattering angles. This 
is the case for the corrected optical depth in Figure 8. Here again a high optical depth 
is noted on the left periphery of the figure. It extends eastward into the upper and cen- 
tral areas of the figure. This may be due to dust intrusions from Asia. Another possi- 
bility for the increasing optical depth to the west is related to the retrieval process. This 
is the result of the satellite sensor's aperture moving to the left and to the right of the 
of the satellite's path to view side regions. The satellite is looking mainly at backscatter 
when it looks directly at the Earth. As the sensor looks to the side it is no longer viewing 
backscatter but smaller scattering angles. This results in a decrease in the phase function 
and a corresponding increase in the optical depth. This is an area that should be inves- 
tigated further in future work. 

Figure 9 shows Benedict's (1989) composite optical depth averaged geographically 
for all longitudes. There is a peak in the vicinity of40°N indicating a potentially greater 
number of particles in the region. This could be a reflection of the advection of Asian 
particulates into the area. These could be natural dust or anthropogenic pollutants. 
There is another peak near 8°X. Since Figure 8 revealed little optical depth from the 
eruption of Kilauca, this peak is probably due to DMS production in the tropics. 

The optical depths obtained from this study, shown together in Figure 10, reveal a 
profile with greater structure, especially in lower latitudes where there is a peak at 
22°X not revealed by Benedict (1989). The peak at 6°X is also shifted north to near 10 
°X. This shows possible error in the 1 1G phase function parameterization. The increase 
in optical depth south of 20°S is subject to error as it is near the data set boundary. The 
values of optical depth are larger than those of Benedict due to the ozone correction. 
This correction increases channel one radiance resulting in a higher optical depth. 

Figures 11 and 12 present the composite particle size index and particle size index 
corrected for water vapor for the entire experimental period. The uncorrected particle 
size index shows a plume in red extending westward from the Kilauea volcano. The 
higher particle size index indicates smaller particles. To the north in the midlatitudes is 
a darker area of blue indicating larger particles. This could be attributed to dust from 
the Asian continent being advected to the east. The large area of dark blue to the 
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Figure 7. Composite channel one optical depth: the area extends from roughlv 50 

°N to 1 2°S and 162"B to I42°W. 
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Figure 8. Dry composite channel one optical depth: the area extends Horn rouehlv 

50° N to 12 ° S and 162°E to 142°W. 
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Figure 9. Composite channel one optical depth (Benedict 1989) 
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Figure 10. Composite uncorrected and dry channel one optical depth: the uncor- 

rected values are smaller than the dry values. 
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western side of the image is unidentified although it could represent more dust. The 
somewhat lower particle size index on the left side of the figure might also be related to 
the retrieval process discussed above. 

The particle size index corrected for water vapor, Figure 12, shows lower values. 
Since water vapor causes an overestimation of S l2 , this was expected. Now the plume 
from Kilauea is less evident but the dust influence to the north is more significant. 

Figure 13 shows the composite particle size index obtained by Benedict (1989). The 
results are averaged geographically for all longitudes. A peak in the neighborhood of 
30-40°N is evident indicating possible influence by smaller particles being advected from 
the west. Another peak is noticed at about 6°N indicating smaller particles, possibly a 
result of DMS production in the tropics. 

The plots of composite particle size index and dry particle size index for this study 
are both shown in Figure 14 for comparison. They show more structure and definition 
than the plot of Benedict (1989). The values are higher than those of Benedict (19S9) 
due to the modifications in the algorithm. In Benedict's (1989) study ozone was over 
corrected for leading to a proportionately higher radiance in channel two and thus lower 
particle size indices. The uncorrected index is higher than the corrected index as before. 
Both curves show the midlatitude peak of Benedict (1989) but also show more of a peak 
in the tropical regions. The Benedict (19S9) plot decreases toward the equator while the 
plots of this study increase. This supports the changes to the algorithm since the particle 
size index is expected to increase in the tropics because of higher DMS production. In 
addition, the dry index has a larger correction for the tropical regions due to the greater 
effect of water vapor there. The increase in S 12 near 60°N may be due to a reduction of 
valid data near the edge of the area. 

The next data sets are in an area of plus or minus five degrees of longitude from the 
ship's track and within plus or minus three days of the ship's position which limits the 
data used in the north-south direction. This is to get a better comparison with the in- 
situ measurements. 

Optical depth near the ship's track from Benedict (1989) is shown in Figure 15. A 
large double peak at about 40°N is evident with a smaller peak at about 33°N. In Figure 
16 this study show's the same peaks at midlatitudes but also a large peak at 22°N. This 
is most likely due to volcanic particles which may have been averaged out in the previous 
composite (Figure 10) w'hich included the entire data set. The corrected optical depths 
on the same plot are even larger. Higher optical depths were expected due to both the 
water vapor correction to the particle size index and the ozone correction. The optical 
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Figure II. Composite particle size index: the area extends from rouehlv 50° N to 

12°S and 162°E to I42°\V. 
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Figure 13. Composite particle size index (Benedict 1989) 
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Figure 14. Composite uncorreded and dry particle size index: The uncorrected 
values are larger than the corrected values. 
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depth measured by a hand-held sun photometer on the ship is also plotted on Figures 
15 and 16. Note that optical depth is scaled differently in the two figures. The magni- 
tude of the optical depths obtained in this study appear to be supported by the 
photometer measurements. This indicates support for the corrections made to the al- 
gorithm. There appears to be a northward shift of about 5-7° in the satellite estimated 
optical depth compared to the ship measurements. This could be caused in part by the 
way the data is averaged by three days around the ship's track. It is possible that sat- 
ellite data within one day of the ship's track would produce better results. The reason 
this wasn't done was to ensure enough data were available to fill in times when cloudy 
pixels were thrown out. A second reason for the shift may be due to inaccuracies in the 
positioning of either the ship or satellite data and a third alternative is that the 
parameterization of the MG phase function is incorrect. Changes to improve the phase 
function are discussed below. 

In-situ measurements of condensation nuclei are shown in Figure 17 although they 
too are shifted several degrees south of the satellite data. The condensation nuclei are 
not the same as CCN, CCN are a subset of the condensation nuclei, being the 
condensation nuclei capable of nucleating with water vapor to form cloud droplets. 
However, for the purpose of optical depth measurement they can be treated similarly. 
Another peak not resolved by Benedict's (1989) or this data set is evident at 6 and 
12°S. 

Figure 18 from Benedict (1989) shows the particle size index in the neighborhood 
of the ship's position. It is a little more definitive than the previous composite and 
reveals a midlatitude peak around 30-40°N, another at about 20°N, and another near 
10°N. A distinction can be made between an influence from Kilauea volcano at 20 6 N 
and DMS production leading to non-sea-salt sulfate aerosols near the equator. 

Figure 19 shows the combined uncorrected and dry particle size index near the re- 
search vessel's track for this study. The peaks are more pronounced than Benedict 
(1989) and there is a minimum at 22°N. where Benedict (19S9) shows a peak. The values 
of S 12 are higher overall and the difference between uncorrected and dry particle size in- 
dex is larger near the equator. There is also an increase in the results moving south of 
the equator as well as a peak not seen in the Benedict (1989) results, possibly due to 
DMS production. 

Next, the non-sea-salt sulfate data obtained from the R, V Oceanographer are plotted 
in Figure 20. Here there is a peak clearly evident centered at about 20°N. This is strong 
indication of smaller particles, most likely from the volcanic fallout. These results point 
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Figure 15. Channel one optical depth along ship's track (Benedict 
1989): photometer measurements superimposed. 
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Figure 16. Uncorrected and dry channel one optical depth along ship's 
track: photometer measurements superimposed 
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Figure 17. Condensation nuclei along ship's track (Benedict 1989) 
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Figure 18. Particle size index along ship's track (Benedict 1989) 
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Figure 19. Uncorrected and dry particle size index along ship's track: the uncor- 

rected values are larger than the dry values. 
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Figure 20. Non-sea-salt sulfate along ship's track (Benedict FS 0 ) 



to errors in the phase function parameterization since there is no corresponding peak in 
particle size index as would be expected. 

The next plot to be studied is the geometry variable used to examine the correlation 
between optical depth and phase function. By plotting the normalized measured radi- 
ance as a function of scattering angle for various particle size indices, any correlation 
between the phase function and particle size index is separated out. This can be seen 
by rearranging equation 2.4 to give 



P = J^L 

V><Fo 



= P {Q)5. 



The left hand side represents values that are all measured or known while the two terms 
on the right must be parameterized (the phase function) or calculated (the optical 
depth). By using the optical depth that best fits the data for a given measured particle 
size index, the phase function can be adjusted and improved. The phase function used 
by Frost (1988) and Benedict (1989) was simply one way, not necessarily the best way 
for this data, of parameterizing the phase function. 

The results of P versus 0 are shown in Figure 21 and can be compared to Figure 2, 
the FIG phase function. To make an accurate comparison, the I1G phase function used 
by Frost (1988) should be multiplied by a mean optical depth. This is done for two 
values of the particle size index in Figure 22 and allows for a direct comparison between 
the HG phase function and study results. For a particle size index of 1.7, the MG phase 
function was multiplied by a mean optical depth of 0.2 and for a particle size index of 
2.1, the HG phase function was multiplied by a mean optical depth ofO.l. These values 
of mean optical depth are only approximate. Further evaluation of the study results 
would be required for more accurate values. 

The shapes of both curves are similar but in the area of 0-20° scattering angle the 
results of this study are up to three orders of magnitude lower. This is due partly to the 
observations in the 0-60° bin. Most of the pixels in this bin have a scattering angle be- 
tween 50 and 60° resulting in a P weighted toward the lower values of the higher scat- 
tering angles. The total number of pixels evaluated for this bin was also less than the 
other bins potentially contributing to an unrepresentative P in this area. For the particle 
size of 1.7 in this study, there also appears to be a shift in minimum P values toward 
higher scattering angles than for the HG phase function. 
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By adjusting phase function parameters a better fit can be obtained. In Figure 23, 
after multiplying the HG phase function by a mean optical depth of 0.15 and using the 
S 12 = 1.7 data curve, the following values from equation 2.8 were obtained: a was set 
to 0.95, g, and g 2 were set to 0.5. In this way the form of a, g,, and g 2 can be determined. 

Another aspect of the problem is that the results were always obtained using the 
HG phase function parameterization for channel one. Channel two results are also 
available but they are computed using the phase function for channel one. The phase 
function is wavelength dependent and could result in some error for channel two calcu- 
lations. Using the same procedure as above but for channel two, it should be possible 
to derive a channel two correction to the Henyey-Greenstein phase function. This would 
result in knowing all parameters in the simplified RTE except optical depth, which could 
then be calculated directly for either channel. A new particle size ratio of optical depths 
from channels one and two would then more directly measure particle size information. 
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Figuie 21. Normalized radiance (P) versus scattering angle (0): for various S 13 . 
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Figure 22. Scaled FIG phase function and P versus O: for particle size index of 
1.7 (solid line) and 2.1. 
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Figure 23. Scaled and adjusted HG phase function and P versus 0: for a particle 
size index of 1.7. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



This study examined an algorithm that extracts aerosol particle size index and op- 
tical depth from satellite observations of upwelling radiance over clear ocean areas. By 
utilizing the difference in sensitivity between channels, particle size information is 
measured and optical depth can then be calculated. Data retrieved from the NOAA-9 
AVHRR sensor was reanalyzed and compared to shipboard measurements made by the 
R/V Oceanographer during the RITS-8S cruise. In this fashion the satellite retrieval 
technique is verified and our understanding of marine aerosols on a regional scale is in- 
creased. 

This study improved upon the results of Benedict (1989) after making corrections to 
account for variable earth-sun distance on solar irradiance, making an adjustment for a 
previous over correction to ozone absorption and correcting the computer code for a 
mathematical error. These combined corrections had a greater effect on the results than 
anticipated. The direction of the corrections was positive, however. An improvement 
to the algorithm, which takes into account the absorption of column water vapor, was 
examined and also found to improve the results. 

Finally, the correlation between aerosol optical depth and the Henyey-Greenstein 
phase function was examined. It was concluded from the plots of P versus 0 for various 
particle size indices that some correlation did exist. An adjustment to one case of the 
HG phase function was made with excellent results. It is recommended to calculate a 
new set of phase functions that fit the experimental data. This would account for the 
correlation between optical depth and scattering angle thus making optical depth and 
phase function independent of one another. Then optical depth could easily be calcu- 
lated and a ratio of optical depths would be a new, more accurate particle size index. 
This is because the slope of the particle size distribution, which is part of optical depth, 
would be more directly compared. 

Another recommendation of this study is to separately calculate a correction to the 
FIG phase function for channel two. Then the channel two phase function could be used 
in the calculations of channel two optical depth instead of the channel one phase func- 
tion. It is also recommended that the simplified RTE be adjusted for multiple scatter. 
The single scattering approximation is good up to an optical depth of about 0. 1. The 
results obtained show that the optical depths were higher than 0.1 leading to errors in 
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the results. Other regional studies should be undertaken to continue examining and 
improving the algorithm and its performance. Global scale studies are also useful but 
it is less likely that there will be in-situ data available for verification. 

While this procedure appeared effective in revealing variations in aerosols in the re- 
gion of interest for known sources, it appeared to have only moderate sensitivity to the 
NSS-SO^'s of interest. This makes studying DMS and its role in global climate regu- 
lation more difficult. It is therefore important that improvements to the algorithm 
continue to be made. As discussed by Benedict (1989), studying areas where there is less 
interference from natural or human-made sources is one possibility as is studying regions 
of higher phytoplankton activity, such as the tropics. Additional satellite-inferred aer- 
osol studies are required to guide future algorithm improvements and to provide crit- 
ically needed measurements of global aerosol behavior. 
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